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Abstract 

We study numerically the existence and character of bound states for positive and negative point 
charges shielded by the response of a two-dimensional homogeneous electron gas. The problem is 
related to many physical situations and has recently arisen in experiments for impurities on metal 
surfaces with Schockley surface states. Mathematical theorems ascertain a bound state for two- 
dimensional circularly symmetric potentials V{r) with drrV{r) < 0. We find that a shielded 
potential with drrV{r) > may also sustain a bound state. Moreover, on the same footing we 
study the electron-electron interactions in the two-dimensional electron gas finding a bound state 
with an energy minimum for a certain electron gas density. 

PACS numbers: 73.20.Hb, 71.15.Mb 
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I. INTRODUCTION AND MOTIVATION 



The properties of the two-dimensional (2D) electron gas and in particular the phenomena 
induced by isolated or clustered impurities embedded in it have attracted a large volume of 
experimental and theoretical research. The interest stems from the increasingly more ideal 
realisations of the 2D electron gas, for example, at interfaces of semiconductor heterostruc- 
tures, at semiconductor surface inversion layers, at (noble) metal surfaces as Schockley 
surface states or as quantum well states in metallic overlayers on insulators or on other 
metals. At the same time the development of different photoelectron and scanning tun- 
neling spectroscopies has enabled an increasingly more accurate characterization of these 
systemsii^. 

Many experiments deal directly or indirectly with the existence of bound electron states in 
systems interacting with the 2D electron gas. The impurity-induced electron localization has 
been studied in the Schockley surface state systems^i^. The quantum diffusion of hydrogen 
on metal surfaces reflects the coupling of the hydrogen with the metallic (possibly 2D) 
electron gas^. Finally, a bound state between two electrons in a 2D electron gas has been 
proposed as an alternative pairing mechanism to the phonon coupling for high-temperature 
superconductivity^"^. Remarkably, it was explicitly pointed out^ that the many-body ground 
state of a dilute gas of fermions is unstable to pairing if and only if a two-body bound state 
exists. 

The background for the above-mentioned physics is laid by basic mathematical theorems. 
First, in two dimensions, for any everywhere attractive circularly symmetric potential, V{r), 
there exists always a bound state no matter how weak the potential is. For the case of a 
shallow potential valley characterized by the drrV{r) < condition, an explicit expres- 
sion for the binding energy, E^, was derived in the textbook by Landai>i^. It turns out 
that the binding energy depends exponentially on the inverse of a negative constant (given 
by the condition) in the investigated t(;eaA;-coupling limit. The theorem of Simon^^ extends 
the conditions for a bound state to the drrV{r) = case, i.e., to the not everywhere 
nonpositive or not everywhere nonnegative suitable potentials. 

In this work we present a numerical analysis for the existence of a 2D bound state 
in potentials fulfilling the J drrV{r) = condition. The physically motivated effective 
potentials, V"e//(r) = AV{r), will refer to the perfectly shielded fields of embedded attractive 
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or repulsive unit-charges in a 2D electron gas, and A plays the role of a convenient coupling 
constant to detailed numerics. In the repulsive case the screening is constrained by the fact 
that the maximum of the surrounding hole density is the uniform electron gas density. Beside 
the cases based on the mentioned standard condition, the overscreening and underscreening 
of unit-charges (with A = 1) shall be investigated as well. The corresponding potentials could 
mimic the shielding-dynamics in, for example, standing wave generation^ on surfaces in the 
presence of impurities. Slightly surprisingly we find a bound state also for the overscreened 
attractive and underscreened repulsive potentials. In order to treat also the electron-electron 
effective interaction we use in the Schrodinger equation the reduced mass /i = 1/2. The 
energy of the ensuing bound electron pair has a minimum at a certain electron gas density. 

The rest of the paper is organized as follows. In the next section. Sec. II, we shall deduce 
our physically motivated model potentials to numerical, 2D bound-state calculations. The 
results obtained are presented on illustrative figures, by considering relevant parameter- 
ranges in coupling, screening-tuning, and density of the electron gas. Finally, Sec. Ill is 
devoted to a short summary. Hartree atomic units, h = = me = 1, will be used in the 
equations and discussion below. 

II. MODELS AND RESULTS 
A. Shielded potentials 

We begin with few mathematical expressions. The 2D Fourier-Hankel (F-H) transforma- 
tion of a F{r) function is 



where Jo{x) is the zeroth-order Bessel function. The inverse F-H transform has the form 



These equations will be used below, in model-potential constructions. 

The field of an embedded charged particle is shielded in the 2D electron gas. Thus instead 




(1) 




(2) 
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of the bare Coulomb form, Vc{q) = '^li'/q, one can write in momentum space 



Viq) = ±vM [1 - An(g)], (3) 

for the shielded field around unit-charges of different signs. Here An{q) is the screening 
density in momentum-space. With a unit-norm An{r) one obtains, via the above Eqs.(l) 
and (3), the V{q ^ 0) oc g limiting behaviour if the real-space density decays faster than 
r""^. This case corresponds to the J drrV{r) = condition of Simon's theorem. 

Notice, that the well-known^^ quasiclassical, Thomas-Fermi (TF) approximation for im- 
purity screening in a 2D electron gas results in a monotonic potential 

Negative and positive signs in Eq.(4) refer to embedded attractive and repulsive unit charges, 
respectively. At large distances (r oo) the Vrrir) potential falls of as ±l/(4r^), and one 
gets the drrVrpir) = ±(1/2) condition. Therefore, this type of potential with negative 
sign belongs to the class analyzed by LandaiJ^. 

We shall characterize the An{r) density directly by the properly normalized hydrogenic 
(H) and gaussian (G) forms, (a^/27r) exp(— ar) and (/^^/vr) exp(— /?^r^), respectively. The 
corresponding Fourier- Hankel transforms are Q;^/(a^ + q^Y^"^ and exp[— g^/(2/5)^], respec- 
tively. Using the former in Eq.(3) and applying Eq.(2), we get 

Vuir) =±^-{l-2u' [h{u)K^{u) - h{u)K,{u)]) , (5) 
in which u = ar/2 for shorthand. The gaussian-form results in 

Voir) = ±^(l-v^/o(^)e-^), (6) 

where z = /5^r^/2. In the above equations Ii{x) and Ki{x) are modified Bessel functions. 

We stress that V{q = 0) = 0, i.e., drrV{r) = in these models. For long distances 
these shielded potentials decay as =F3/(2a^r^), and =Fl/(4/?^r^), respectively. The attractive 
potentials are exhibited in Fig. 1, by fixing a = 4 (solid curve) and p = 2^2 (dashed curve . 
The a = 4 value could correspond to a strictly atomistic, Is-like in 2D, electron density. 
For a detailed comparison, two other potentials are also plotted. That is, the bare Coulomb 
one, Vc{r) = (— 1/r), and the conventional Thomas-Fermi form, Vrpi^) of Eq-(4). 
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FIG. 1: (Color online) The shielded attractive potentials of Eq.(5) (solid blue curve) and Eq.(6) 
(dashed red curve) with a = 4 and /3 = 2\/2, respectively. The dotted black curve refers to a bare 
Coulomb potential, — 1/r. The dash-dotted green curve is devoted to the attractive Thomas-Fermi 
form of Eq.(4). 

In order to go beyond the drrV{r) = condition, one may multiply the shielding 
parts of Eqs.(5) and (6) by a A variable [See Eqs. (fTUj) and (fT^ below]. In such a way 
we can investigate the overscreened (A > 1) and underscreened (A < 1) cases, for which 
Jq°° drrV{r, A) 7^ 0. Moreover, as discussed above the effective potentials have a multiplica- 
tive coupling constant A, [Ve//(r) = AV{r,X)]. In our 2D numerics A and A will serve as 
convenient parameters. 




B. Numerical solution of the 2D Schrodinger equation 

The bound-state energy levels (Ef,) and wave functions ['?/'(r)] satisfy the 2D Schrodinger 
equation 



m = 0, 



(7) 



where /i is the reduced mass; it is unity for the impurity-electron case. 
In circular symmetry the wave function separates as 



^(r) 



2tt 



Rmnir) 



(8) 



where m = 0, ±1, ±2.... is the azimuthal quantum number and n = 1, 2, 3, .... is the radial 
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quantum number related to the number of radial nodes (n — 1) of the radial wave function 
Rmn{f)- In this work we are only interested in the values m = and n = 1. Further, by 
making the substitution Umnir) = r^/'^Rmn{f) we obtain the differential equation 



This is the same form as the radial equation studied in spherically symmetric problems. 
We solve the equation on an exponentially expanding radial mesh, r(j) = Tmin exp((j — 1) Ax) 
with j=l . . . N. With a given guess for the eigenvalue Ef, the function Umn{f^) is integrated 
outwards from origin and inwards from a large radius by starting with its asymptotic expan- 
sions. At a matching point close to the classical turning point the logarithmic derivatives 
of the outward and inward integrated solutions are required to coincide by adjusting the 
eigenvalue Eh- The parameters of the radial mesh, Tmin, Ax, and N are varied until the 
numerical convergence of the eigenvalue is obtained. 

C. Attractive electron-ion interaction 

According to our numerical calculations, the attractive Vtf{^) of Eq.(4) gives = 
-0.2853, while the Vff(r) of Eq.(5) results in E^ ' = -0.0296 values for a =4, in Hartree 
units. There is an about order-of-magnitude reduction in binding due to an atomistic screen- 
ing, beyond quasiclassics. 

In order to study the so-called not everywhere nonpositive case (physically: the shielded 
positive charge) in more detail we will use, without loss of generality, a = 4 in Eq.(5) and 
modify it as 



In Fig. 2 we have plotted the E^ (A) energies obtained with A=l and A G [0.97, 1.03]. The 
A-tuning refers to small under- and over screening. The figure clearly shows the sensitivity 
of binding on the details of shielding. For A > 1 one has drrVeff{r) > 0, with Eq.(|10|). 

Additional information on the A-dependence of E^^\A) for Eq. ljlOp and A=l, are given 
in Fig. 3 with A G [1.25,0.75]. In this case drrVeff{r) = for any (finite) A, as we 
pointed out earlier. In harmony with Simon's theoren>ii, |Sb(A)| < exp(— cA^^), we get for 
our case [Ej^^\A)] an about c = 5.292 value for the suitable constant in the investigated 





Vfffir) = -A-(l-8Ar2[/o(2r)Ki(2r)-/i(2r)iro(2r)]). 



(10) 
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FIG. 2: Binding energy, '{X), based on Eq.(7) with Eq.^ for A=l and A E [0.97, 1.03]. 

coupling-parameter range. The macroscopic dielectric constant of a real medium can result 
in reduced, A < 1 values to effective interactions. 

In this attractive impurity case we have shown a remarkable sensitivity of theoretical 
bound-state characteristics to shielding conditions. This is in accord with important spec- 
troscopic information obtained^' by scanning tunneling spectroscopy for different adatoms in 
generated standing wave patterns. The observed peak-shift and amplitude-decrease in dif- 
ferential tunneling conductance, as adatoms are approached to a step on the surface, signal 
the experimental sensitivity. 

We note that our detailed numerical analysis is based on an effective Schrodinger equation. 
Further attempts are needed therefore to consider the many-body aspects of the localization 
problem in more detail. For example, the proper description of the width of an adatom- 
induced weakly bound state and its influence on scattering characteristics are important 
questions to future studies. 

D. Repulsive electron-electron interaction 

The so-called not everywhere nonnegative case (physically: the shielded negative charge) 
will be treated similarly via Eq.(5) of positive sign, but with the constraint 



in order to model bounded, complete depletion of the 2D electron gas density at r = 0. Here 
is the density-parameter of the 2D gas with density no- The same constraint is used to 



AaV27r = tiq = l/{nrl) 



(11) 
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FIG. 3: Binding energy, Elf^\A), based on Eq.(7) with Eq.(lO) for A=l and A G [1.25,0.75]. 



the corresponding (positive sign) reparametrized Eq.(6) with 



A/^Vvr = no = l/{nrl) 



(12) 



as 



+ (l - XV2^ Io{z) e-'^ , (13) 

where, according to the above discussion, z = (r/rs)^/(2A) now. Remember, that at A = 1 
one has /q°° drrVeff{r) = 0, independently of the value of a finite A. 

As we use the depletion-constraint-based densities to model effective electron-electron in- 
teractions, it is important to give an additional physical argument on their proper behaviour. 
The interaction energy (e) of a repulsive point charge with the surrounding, normalized 
(A = 1) shielding hole can be characterized by the classical equation 



1 f°° 1 

e = 271 drr-An{r). 

2 Jo 



(14) 



This equation results in the €h = — l/(v^'"s), and sg = —V^K'^'^s) expressions for the 
exponential and gaussian screening, respectively. These energies are between the values 
based on the exchange-only, = — 4v^/(37rrs), and Wigner-monatom ,^^i^^ ew = —^/^s, 
limiting approximations. The latter corresponds to the An{r < r^) = Uq extremum model 
for the hole-density, while the former to the Pauli-hole^^ of an ideal 2D system. 

The two repulsive potentials, discussed above in details, are used in Eq.((7j) with the 
reduced mass fi = 1/2, and the following numerical results are obtained. The binding 
energies Eb{rs) are plotted in Fig. 4, as a function of the density parameter r^. Computed 
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FIG. 4: (Color online) Binding energies, Ei,{rs), in the shielded field of a negative charge. Data 
are based on Eq.Q with /x = 1/2, and A = A = 1 in the applied two effective interactions. The 
solid blue and dashed red curves refer to hydrogenic [Eq. (5)] and gaussian [Eq. (6)] models, 
respectively. See the text for further details. 

values are connected by solid (the case of exponential shielding) and dashed (the case of 
gaussian model) curves. The other parameters behind these results are fixed as A = A = 1. 

Remarkably, there are Vg parameter values at which the binding energies are optimal, i.e., 
they have extremal values. By increasing or decreasing the density of the electron gas, the 
bindings become weaker. The extremal values are at about = 11.4 for the hydrogen-like 
model, and at = 17.2 for the gaussian model. The binding energies are in the 10~^ range 
in atomic units, for the dilute system. The low carrier density and thus a small Fermi energy, 
Ep = l/'^s, are important characteristics of cuprate superconductors^^. In these materials 
the ratio of the critical temperature (Tc) and the Fermi energy is in the range of 10~^. 
Furthermore, there is a saturation and supression of Tc with increasing carrier density; for 
further detail we refer to Fig. 3 of Ref.[16]. 

The observed extremal character is related, physically, to our pseudononlinear construc- 
tion of effective shielding of repulsive charges. Namely, the constraint via the bounded 
depletion hole at r = fixes a reasonable scaling in the potentials. In standard linear- 
response theory^ the hole-density at r = can become higher than the host density no- 
In such an attempt one can get a monotonic dependence of Ef,{rs) on in the resulting 
shielded fields. 

Further information is given in Figs. 5 and 6, which show the 2D radial densities. 
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FIG. 5: The radial density 27rr|V'(?')P and the potential based on hydrogenic screening of a negative 
unit-charge. These are computed at the = 11.4 value of the density parameter. 
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FIG. 6: The radial density 2^^r\^p(r)\'^ and the potential based on gaussian screening of a negative 
unit-charge. These are computed at the = 17.2 value of the density parameter. 

27rr|-?/'(r)p, computed with bound-state wave functions at the extremal Wigner-Seitz pa- 
rameters. The corresponding potentials are also plotted. As expected, the square-integrable 
wave functions are localized at about the potential minima. In the investigated equal-mass 
case this extension can represent a certain coherence- lenght; somewhat surprisingly it is only 
twice of the extremal density parameter. 

As we observed in Fig. 3 for the attractive case, the binding energy depends on the 
coupling, A. An illustration on this fact for the repulsive case is given in Fig. 7. The 
gaussian model of Eq. ()13|l is used with A = 1, /i = 1/2, and = 17.2. Results for A < 1 
are plotted. We can approximate our data by a quadratic expression, i.e., \eI'^\A)\ ~ A^. 
Rescaling of A = 1 by a macroscopic dielectric constant could result in a notable reduction 
of the above-mentioned (see. Fig. 4, also) binding energies. 

We finish our representation, by discussing the question of undershielding the repulsive 
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FIG. 7: The coupling-constant (A) dependence of the binding energy for the repulsive case. The 
results are based on the gaussian model of Eq. (|TH|) with /x = 1/2 in Eq.Q. The Wigner-Seitz 



parameter is = 17.2 and A = 1. 
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FIG. 8: The shielding-constant (A) dependence of the binding energy for the repulsive case. The 
results are based on the gaussian model of Ea.lfT^ with fi = 1/2 in Eq.Q. The Wigner-Seitz 
parameter is = 17.2 and A = 1. 

charge. Fig. 8 is devoted to this problem. We have used Eq. (fT^ with A = 1 and A < 1 in the 
Schrodinger equation, Eq.Q with fi = 1/2, at = 17.2. One can observe (see. Fig. 2, for 
the attractive case) essential reductions of energies for the drrV{r) > unconventional 
condition. 

The last results of Figs. 7 and 8, together with Figs. 4-6, signal a nontrivial sensitivity of 
the magnitude of the binding energy on the concrete physical situation. Fortunately, there 
is a physical limitation. It should be clear from the discussion at Eq. (fn|l . and as our illus- 
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trative figures indeed show, that a more localized real-space character of the bounded (and 
normalized) depletion hole around a repulsive unit-charge, results in an effective potential 
with a repulsive part of a shorter range. Clearly, the limitation is given by the Wigner 
model, in which one has the An(r < r^) = Uq extremum for the hole. 

III. SUMMARY 

In this work we have investigated the problem of bound states in two-dimensional shielded 
potentials. Effective potentials, based on direct approximation for the screening charge den- 
sities around attractive and repulsive unit charges, are employed. Particularly, the effective 
electron-electron interaction is modelled via a properly constrained depletion hole. In the 
detailed numerical analysis performed, we found that in both of the basically attractive, and 
basically repulsive potentials bound states appear under the drrV{r) = standard, and 
may appear under the drrV{r) > unconventional conditions. In the repulsive case an 
extremal character of the binding energy, as a function of the density of the host 2D electron 
gas, is estabhshed with our physical models for effective electron-electron interactions. 
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